When Did Cosmic Acceleration Start ? 
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A precise determination, and comparison, of the epoch of the onset of cosmic acceleration, at 
redshift z acc , and of dark energy domination, at z eq , provides an interesting measure with which to 
parameterize dark energy models. By combining several cosmological datasets we place constraints 
on the redshift and age of cosmological acceleration. For a ACDM model, we find the constraint 
Zacc = 0.76 ± 0.10 at 95% c.L, occurring 6.7 ± 0.4 Gyrs ago. Allowing a constant equation of 
state but different from —1 changes the constraint to z acc = 0.81 ± 0.12 (6.9 ± 0.5 Gyrs ago) and 
Zcq = 0.48 ± 0.14 (4.9 ± 0.9 Gyrs ago) , while dynamical models markedly increase the error on the 
constraint with z acc = 0.81 ± 0.30 (6.8 ± 1.4 Gyrs ago) and z eq = 0.44 ± 0.20 (4.5 ± 1.0 Gyrs ago). 
Unified dark energy models as Silent Quartessence yield: z acc — 0.80 ± 0.16 (6.8 ± 0.6 Gyrs ago). 



I. INTRODUCTION 

The existence of a dark and unclustered energy 
component responsible for more than 70% of the over- 
all density of our universe is now suggested at high 
significance by most of the latest cosmological data 
(see e.g. Q, 1). 

A cosmological constant provides a possible candi- 
date for the dark energy component, but it needs to 
have its initial conditions properly 'tuned' in order 
to dominate the universe expansion at precisely the 
present time. Indeed the energy density in a cosmo- 
logical constant pa does not evolve, while both matter 
(p m ) and radiation (/?,.) energy densities evolve rapidly 
with the expansion of the universe. The small current 
value of pa implies an extreme fine-tuning of initial 
conditions with p\/p r — 10" 123 at the Planck time 
(when the temperature of the universe was T ~ 10 19 
GeV), or p\/p r ~ 10~ 55 at the time of the electroweak 
phase transition (T ~ 100 GeV). Moreover, pa lies 
in a very small window, since a slightly larger value 
makes the universe accelerate much before the present 
epoch, thereby inhibiting structure formation, while a 
negative value may cause the universe to re-collapse. 

The cosmological constant problem has, subse- 
quently, motivated several "dynamical" alternatives 
(see e.g. [|| for a recent and very complete review) as 
a slowly-rolling scalar field, "quintessence" or 
a "k-essence" scalar field with non-canonical kinetic 
terms in the Lagrangian Q, string-inspired models 
such as the contribution of nonlinear short distance 
physics to vacuum energy 0] , and modified Friedman 
equations at late time [8J or large distances @. 

Other possibilities include anthropic arguments 
( [ToL [Til [l2 |) and "backreaction" of non linear inho- 
mogeneities (see [l3j], but see also [lil. Il5|). 

It is plausible that a solution to the dark energy 
problem could be found by identifying a time corre- 
lation between the epoch of appearance of this ex- 
otic component and a well understood and physically 
motivated event such as the time of matter-radiation 
equality, the origin of non linear structures or, ulti- 
mately, life. It is, therefore, clear that a first crucial 
measurement that has to be made is the determina- 



tion of the redshift and time of dark energy domina- 
tion. Evidence for dark energy at very high redshifts 
(z > 1), when the cosmological constant is negligi- 
ble, would indeed favor models based on scalar fields, 
possibly coupled to dark matter [l6|, [llj • While the 
appearance of dark energy at lower redshifts (z < 0.2) 
would, on the contrary hint at "phantom" (w < — 1) 
models. Anthropic principle arguments are definitely 
less appealing if dark energy dominates well after 
the epoch of formation of terrestrial planets. At the 
same time, backreaction models could be perceived as 
much less convincing if dark energy starts in a time 
when nonlinear structures are already well present and 
formed. 

The starting point of cosmic acceleration, however, 
is not a model independent quantity. If the universe 
is in accelerated expansion today we can identify two 
crucial epochs. Firstly, the epoch of equality between 
matter and the dark energy component, at redshift 



, defined as 



Pm(Zeq) = Px(z eq ) 



(1) 



where p m (z) and px(z) are the energy densities of the 
matter and dark energy components at redshift z re- 
spectively. This epoch is generally different from, and 
follows, the redshift z acc when the universe started to 
accelerate, defined as 
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where H — a/ a is the Hubble parameter at redshift z. 
The two epochs are model dependent and distinct. In 
the simple case of a cosmological constant, for exam- 
ple, the two redshifts are not equal and the following 
simple relation: 



2 1 / 3 (l + z eg )-l 



(3) 



holds. The age of the universe at each of those red- 
shifts can then be easily computed from: 



t(z) 



dz 



(4) 



(l + z)H(z) 

and compared with the current age of the universe to 



2 



It is clear that constraints on the model-dependent tad and Mortsell (HM), see [24(), where: 



quantities z acc , z eq , t(z eq ) and t(z acc ) can provide rel- 
evant information for several studies. In this paper, 
we focus on constraining these quantities with cur- 
rent cosmological data with the goal of clarifying the 
following points: how model independent are the con- 
straints on the epoch of dark energy domination ? how 
does a different choice of cosmological datasets or pa- 
rameters affect those constraints ? Finally, are the 
constraints, derived in a general dark energy scenario, 
consistent with the predictions of a cosmological con- 
stant ? 

Our paper is organized as follows, in the next sec- 
tion we introduce our data analysis method, describ- 
ing the datasets and dark energy parameterizations 
adopted. In section III we present the results of our 
analysis and in section IV we derive our conclusions. 



II. LIKELIHOOD ANALYSIS 

The method we adopt is based on the publicly avail- 
able Markov Chain Monte Carlo package cosmomc 
. We sample the following dimensional set of cos- 
mological parameters, adopting flat priors on them: 
the physical baryon and CDM densities, ojf, — il^h 2 
and uj c = £l c h 2 , the ratio of the sound horizon to 
the angular diameter distance at decoupling, 9 S , the 
scalar spectral index, n s , and the optical depth to 
reionization, r. Furthermore, we consider purely adi- 
abatic initial conditions and the possibility of curved 
universes, Vt to t ^ 1- We also consider the possibility 
of having a running of the spectral index dn s /dlnk at 
k = 0.002h~ 1 Mpc (see e.g. [131), an extra-background 
of relativistic particles (parametrized with an effective 
number of neutrino species N"^ ^ 3, see e.g. [20]) 
and a non-zero, degenerate, neutrino mass of energy 
density (see e.g. |21|): 



92.5eU 



(5) 



in order to establish how robust measurements of z acc 
and z eq are to broader cosmological models. 
Finally, we will also investigate the possibility of a 
dark energy equation of state w different from —1. 
Other than a constant equation of state w we consider 
the possibility of a varying with redshift equation of 
state. In particular we consider a linear dep endence 
on scale factor a = (1 + z) _1 as [22[ and |23j |: 



w(a) = wq +wi(1 — a) 



(0) 



where the equation of state changes from wq to wq+wi 
at higher redshifts. We refer to this as Chevallicr- 
Polarski-Lindcr (CPL) parameterization. 

We also consider a more sophisticated parametriza- 
tion that takes in to account the rate and redshift of 
the transition. We use the model proposed by Hannes- 



w(a) = wqWi 



w\a q + wqcis 



(7) 



In this model the equation of state changes from Wq 
to w\ around redshift z s = 1 — l/a s with a gradient 
transition given by q. We assume u>o,i > —3, 0.1 < 
a s < 1.0 and 1 < q < 10 as external priors for this 
model. 

Finally, we consider the Quartessence (or Chaply- 
gin gas) model (see (25[) as unified dark energy-dark 
matter model. In this scenario cold dark matter and 
dark energy are the same fluid with equation of state: 



w(a) 



Wq 



~w + (1 + wo)a~ 3(a+1) 



(8) 



where wq is the current value of the equation of state 
and a is a parameter that has to be constrained from 
observations. From the equation above, it is clear that 
at early times, when a — > 0, we have w — > 0, and the 
fluid behaves as non relativistic matter. At late times, 
when o > 1, we obtain w — ► — 1. The matter cluster- 
ing presents strong instabilities and oscillations in this 
model [27| unless one assumes intrinsic non-adiabatic 
perturbations such that the effective sound speed van- 
ishes 26] . In this paper we therefore consider only this 
"Silent" Quartessence. 

The MCMC convergence diagnostics are done on 
7 chains applying the Gelman and Rubin "variance 
of chain mean" /"mean of chain variances" R statis- 
tic for each parameter. Our 1 — D and 2 — D con- 
straints are obtained after marginalization over the re- 
maining "nuisance" parameters, again using the pro- 
grams included in the cosmomc package. Tempera- 
ture, cross polarization and polarization CMB fluctu- 
ations from the WMAP 3 year data 0, [H [H, HI 
are considered and we include a top-hat age prior 
10 Gyr < t Q < 20 Gyr. We combine the WMAP data 
with the the real-space power spectrum of galaxies 
from the Sloan Digital Sky Survey (SDSS) [3l| and 
2dF survey We restrict the analysis to a range 

of scales over which the fluctuations are assumed to 
be in the linear regime (technically, k < 0.2/1" 1 Mpc) 
and we marginalize over a bias b considered as an ad- 
ditional nuisance parameter. We also incorporate the 
constraints obtained from the supernova (SN-Ia) lu- 
minosity measurements by using the so-called GOLD 
data set from [32| and the Supernovae Legacy Survey 
(SNLS) data from [H. 



III. RESULTS 

A. Cosmological Datasets 

Using the analysis method described in the previous 
section we have constrained the value of z eq , t eq , z acc 
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Dataset 
WMAP+ 


Zeq 


to — teq 

[Gyrs] 




to t a cc 

[Grys] 


to 

[Gyrs] 


Alone 

+SDSS 

+2dF 

+GOLD 

+SNLS 

+all 


n 47+0-09 

u '^' -0.09 

n 4n+ 008 

SS+ 06 

u '°°-0.06 
U - 4O -0.06 

n 4o +0 - 04 


4.7l° j 
4 s+ 05 

4 -8- .3 
4 1+0-4 

4 6+°- 4 

^■"-0.4 

4 s+o- 3 


0.861°- 
0.771°;!° 
87+ ' 07 

n 74+0.08 
U - ' 4 -0.08 

n QQ+o-08 
76+ 05 

u -'°-0.05 


7 O+O' 4 

6 7+0- 3 

D -'-0.3 

7 1+0.2 
'• L -0.3 

6-6±g:l 

6.9l° j 


13.8«- 3 3 
13.8±g-i 
13.8l°,; 2 
13.S±°,; 2 2 

13.81°;! 
13.9l°;J 



TABLE I: Constraints on z eq , t eq , z acc and t aC c, at 68% 
c.l., in comparison with various datasets for ACDM. 
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to — teq 

[Gyrs] 


^acc 


to t a cc 

[Grys] 


to 

[Gyrs] 


Q.tot + 1 

dn/dlnk =fi 

K tt + 3 
Sm„ > 


0.48l°;° 7 
S2+ 010 

u - oz -0.10 
n oy+0.05 
u -°' -0.05 

40+ 05 

u -^ u -0.06 

S7+ 04 

U - J ' -0.04 


o n+0.8 
J - M -0.8 
4 1+0-3 

4 s+o- 5 
4 n+O-3 

*' Z -0.3 


0.811^6 
0.68l° :i ° 

n 79 +0.06 
u - ' z -0.10 
n 77+O.O6 
u -' '-0.06 
n 70+O.O5 
u -'°-0.05 


6.9l°j 

6 q+°- 3 
6.6l°; 2 2 

6.8±g : | 

6.71°;?, 


13.9l°; 2 

IS 1+0-8 

ro.i_o.9 
14 1+ - 1 

14 l+°' 2 



TABLE II: Constraints on z eq , teq, z acc and tacc, at 68% 
c.l., under differing theoretical assumptions for the under- 
lying cosmological model. 



and t acc in light of the various datasets and cosmolog- 
ical scenarios. This kind of test is extremely useful in 
order to identify the presence of possible systematics. 

The constraints on z eq , t eq , z acc and t acc for various 
datasets are reported in Table I. 

As we can see, there is a general agreement between 
the results: namely, in a cosmological constant model, 
dark energy became the dominant component at red- 
shift z eq ~ 0.4, 4.3 Gyrs ago, and the accelerated ex- 
pansion of the universe started at z acc ~ 0.75, 6.7 Gyrs 
ago. Since we are assuming a cosmological constant 
z eq and z acc are not independent but follow Eq.(3). 
It is interesting, however, to note that the SDSS and 
GOLD datasets seem to favour a lower redshift for the 
dark energy's dominant than that suggested when the 
2dF or SNLS datasets are included, respectively. 
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to — t eq 
[Gyrs] 
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[Gyrs] 


w # -1 
CPL 
HM 
SQ 


4S+ 07 

U - 4J -0.06 

44+ - 11 
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4 S+0- 5 
4 '°-0.5 
4 Ci+0-7 

4 6+°- 6 


0.79l°;°7 
0.80!°;!? 
79+°- 14 

<J- li> _ .14 

n 8n+°- 08 
0-80_o.o8 


6.8l° j 
6.81°;? 

6.7±g:S 
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13.8l°; 2 
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TABLE III: Constraints on z eq , t eq , z acc and t acc , at 68% 
c.l., for different theoretical assumptions about the nature 
of the dark energy component. 



dard scenario, could provide a hint for the presence of 
an extra background of relativistic particles [35| ■ 



C. Dynamical Dark Energy 



B. Theoretical assumptions on the background 
cosmological model. 

Since the results appear stable to the inclu- 
sion/exclusion of the experimental datasets, we now 
consider the full set of cosmological data and study the 
dependence on some of the theoretical assumptions 
on the background cosmological model. We consider 
possible variations from —1 in a constant dark energy 
equation of state, non-flat universes, a running of the 
spectral index of the primordial inflationary pertur- 
bations, massive neutrinos and an extra background 
of relativistic particles. All those constraints are re- 
ported in Table II. As we can see, the results are 
consistent with those reported in Table I. However, 
as expected, the constraints are in general weaker. 
While including a constant dark energy equation of 
state w 7^ — 1 has little effect, considering a universe 
with spatial curvature generally doubles the error bars 
on all the parameters and results in a lower redshift 
and time of dark energy domination. It is interest- 
ing to note that considering an extra background of 
relativistic particles has a strong effect on the age of 
the universe and of dark energy. This raises an inter- 
esting question of whether the recent discovery of the 
APM 08279+5255 quasar at z = 3.91, whose age of 
2 — 3 Gyr can't easily be accommodated in the stan- 



We now study the sensitivity of the epochs of dark 
energy domination and the onset of acceleration to 
differing dark energy models. We compare models to 
WMAP+2dF+SNLS datasets, this should be consid- 
ered a more conservative choice in comparison to the 
"all" dataset described in the previous section. For 
the Silent Quartessence, however, we consider only 
WMAP+SNLS, as a unified dark energy model, in- 
clude no cold dark matter, we omit the redshift and 
time of equivalence with the baryonic component. The 
constraints are reported in Table III. 

Allowing for an equation of state which is varying 
with redshift can strongly affect the constraints, with 
error bars as large as four times those reported in Ta- 
ble I. However, the mean values are generally consis- 
tent with the previous results, i.e. there is no indica- 
tion for deviations from a cosmological constant. In 
this sense, it is useful to plot the constraints on the 
z a cc — z eq plane as we do in Fig.l. A cosmological con- 
stant in this plane generates a line described by Eq. 
3.. As showed by three contour plots, while including 
a dynamical component leaves the possibility of a dif- 
ferent relation between the two redshifts, the case of a 
cosmological constant is always well inside the lcr c.l.. 

In Fig. 2 we plot the 2a constraints on the deceler- 
ation parameter q in function of the redshift for the 
four different dark energy parameterizations. While 
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FIG. 1: Constraints on the z acc — z eq plane for different 
dark energy parameterizations. From top to bottom w 
constant, Chevallier-Polarski-Linder (CPL) and Hannes- 
tad Mortsell (HM) models. Also plotted (dashed line) is 
the cosmological constant case z acc = 2 1 ^ 3 (z eq + 1) — 1 
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there is a large spread in the values, especially when 
more complex parametrizations such as CPL or HM 
are considered, there is a very good agreement, and 
all the models point towards the same acceleration 
redshift value at z ~ 0.8. 

IV. CONCLUSIONS 



In this brief paper we have presented several con- 
straints on the epoch and redshift of dark energy 
domination and of cosmic acceleration. We have de- 
rived those constraints using different datasets, differ- 
ent theoretical assumptions and different dark energy 
parametrizations. We have found that a redshift and 
epoch of acceleration at z acc = 0.78 and to — t acc — 6.9 
Gyrs and a redshift and epoch of dark energy domi- 
nation start at z acc = 0.43 and to — t acc = 4.4 Gyrs, as 
expected for a flat universe with = 0.7 is in agree- 
ment with all the possible cases considered. More- 
over, despite the large set of models and data ana- 
lyzed, there a very little spread in the best-fit val- 
ues. Curvature, running of the spectral-index, massive 
neutrinos and an extra-background of relativistic par- 
ticles are non-standard cosmological parameters that 
can strongly enlarge the error bars on z eq and z acc 
in case of a cosmological constant model. Allowing a 
constant or dynamical dark energy equation of state 
different from —1 produces similar results, however 
the best fit values are, again, not significantly altered. 
A tension in the derived best fit values appears when 
considering galaxy clustering data from SLOAN and 
2dF and supernovae type la from Riess et al. and 
SNLS datasets separately. However the significance 
of the discrepancy is well below the 2a c.L. As a final 
remark, we like to stress that the analysis presented 
here relies nonetheless in the assumption of a class 
of scenarios. It may be possible to construct more 
complicated cosmological and/or dark energy models 
that could evade the constraints presented here. Fu- 
ture data may well falsify those possibilities as well 
the simple case of a cosmological constant by testing 
the z„ 



z eq relation of Eq. 3. 



FIG. 2: Constraints at 95% c.L on the deceleration param- 
eter q(z) in function of redshift for different dark energy 
parametrizations (see text). 



Acknowledgments We are pleased to thank 
Rachel Bean for help, discussions and useful comments 
on the manuscript. 



[1] D. N. Spergel et al. , |arXiv:astro-ph/0603449| 

[2] U. Seljak, A. Slo sar and P. McDonald, 

|arXiv:astro-ph/0604335l 
[3] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. 

Mod. Phys. D 15 (2006) 1753 [arXiv:hep-th/0603"057] . 
[4] C. Wetterich, Nucl. Phys. B 302, 668 (1988). 
[5] R. R. Caldwell, R. Dave an d P. J. Steinhardt, Phys . 

Rev. Lett. 80, 1582 (1998) [arXiv:astro-ph/970 8069 . 
[6] C. Armendariz-Picon, T. Damour and 



V. M ukhano v, Phys. Lett. B 458, 209 (1999) 
[arXivihep-th /9904075] . 
[7] L.Mersini, M. Bastero-Gil, P.Kanti, 

Phys.Rev .D64:043508 (2001), |hei>ph/0101210| 
| hep-ph/010"oT34] ; M.Bastero-Gil, P.Frampton, 
Phys.Rev.D65: 106002,(2002) , 
M.Bastero-Gil and L.Mersini, 
hys.Rev.D67:103519,(2003), |hep-th/020527l1 . 
[8] K.Freese, M.Lewis, Phys.Lett.B540:l-8,(2002), 



.Mersini, 
hep-th/0110167 ; 



[9] 

[10] 

[11] 
[12] 
[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 



[astro-ph/0201229] . 

G.Dvali , G.Gabadadze, M.Porrati, Phys.Lett.B 485, 

(2000) , [hep-th/0005016] ; G.Dvali a nd G.Gabadadze , 
Phys.Rev.D63,065007, (2001), [hep-th/00080"54] ; 
G.Dvali and M.Turner, [astro-ph/0 301510 . 
Barrow, J.D. and Tipler, F.J. (1986) The Anthropic 
Cosmological Principle, New York: Oxford University 
Press. 

Garriga, J. an d Vilenkin, A., PR D 61 083502 (2000). 
Weinberg, S., |astro-ph/0104482| 

E. W. Kolb, S . Matarrese and A. Riotto , New J. Phys. 

8 (2006) 322 [arXiviastro-ph /0506534] . 

E. E. Flanagan, Phys . Rev. D 71, 103521 (2005) 

[arXiv:hep-th/0503202| . 

C. M. Hirata and U. Seljak, Phy s. Rev. D 72 (2005) 
083501 [arXiv:astro-ph /0503582] . 
L. Amendola, Phys. R ev. D 62, 043511 (2000) 
[arXiv:astro-ph/9908023] . 

R. Bea n and J. Magueijo, Phy s. Lett. B 517, 177 

(2001) [arXiv:astro-ph/0007199l 
A. Lewis and S. Bridle, Phys 
103511 [arXiv:astro-ph/0205436 
W. H. Kinney et al., Phys. Rev. 
arXiv:astro-ph/0605338 ; W. H 



Rev. D 66 (2002) 

D 74 (2006) 023502 
Kinney et al, Phys. 



Rev. D 69 (2004) 103516 arXiv:hep-ph/0305130 . 
R. Bowen et al, Mon. Not. Roy. A"stron. 
Soc. 334 (2002) 760 [arXiv:astro-ph/0110636] ; 
S. H. Hansen et al, Phy s. Rev. D 65 (2002) 023511 
[arXiv:astro-ph/0105385l . 
G. L. Fogli et al, Phys. 



[22 
[23 
[24 
[25 
[26 

[27 

[28 

[29 

[30 

[31 
[32 
[33 

[34 
[35 



[arXiv:hep-plT/ 0608060 ; G. L. Fogli et al, P hys. Rev. 
D 70 (2004) 113003 [arXiv:hep-ph/0408045| . 
M. Chevallier and D. Polarski, Int. J. Mod. Phy. D., 
10, 213 (2001) 

E. V. Linder, Phys. Re v. Lett. 90 (2003) 091301 
[arXiv: astro-ph/02085 12] . 

S. Hannestad and E. Mor tsell, JCAP 0409 (2004) 001 
[arXiv:astro-ph/0407259] . 

N. Bilic, G B. T upper and R. D. Viollier, Phys. Lett. 
B 535 (2002) 17 [arXiv:astro-ph/0111325l . 
R. R. R. Reis, I. Waga, M. O. Calvao and 
S. E. Joras, Phys. Re v. D 68, 061302 (2003) 
[arXiv:astrc^ph7 0306004 . 
H. Sandvik, M. Tegmark, 



I. Waga, Phys. Rev. 
[arXiv:astro-ph/0212114] . 
L. Page et al. 
|arXiv:astro-ph/0603450l 
G. Hinshaw et al. 



D 



M. 
69, 



[WMAP 
[WMAP 
[WMAP 



Zaldarriaga and 
123524 (2004) 

Collaboration] , 

Collaboration] , 

Collaboration] , 



Rev. D 75 (2007) 053001 



arXiv:astro-ph/0603451 
N. Jarosik et al. 
|arXiv:astro-ph /0603452] 
M. Tegmark et al., Astrophys. J. 606, 702 (2004). 
A.G. Riess et al, Astrophys. J. 665R, 607 (2004). 
S. Cole et al, Mon. Not. Roy. Astron. Soc. 362, 505 
(2005). 

P. Astier et al, Astron. Astrophys. 447, 31 (2006). 
D. Jain and A. Dev, Ph ys. Lett. B 633 (2006) 436 
[arXiv:astro-ph/0509212] . 



